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A molecular hydrogen absorber at a lookback time of 12.4 billion years, corresponding to 10% of 
the age of the universe today, is analyzed to put a constraint on a varying proton-electron mass 
ratio, p. A high resolution spectrum of the J1443+2724 quasar, which was observed with the Very 
Large Telescope, is used to create an accurate model of 89 Lyman and Werner band transitions 
whose relative frequencies are sensitive to p, yielding a limit on the relative deviation from the 
current laboratory value of A p/p = (—9.5 ± 5.4 ata t ± 5.3 sys ) x 10 -6 


PACS numbers: 98.80.-k, 06.20.Jr, 98.58.-w, 33.20.-t 

The accelerated expansion of the universe is ascribed 
to an elusive form of gravitational repulsive action re¬ 
ferred to as dark energy [I] . Whether it is a cosmological 
constant, inherent to the fabric of space-time, or whether 
it may be ascribed to some dynamical action in the form 
of a scalar field </> [2], is an open issue. In the latter case 
it has been shown that the interaction of the postulated 
quintessence fields <f> to matter cannot be ignored, giving 
rise to a variation of the fundamental coupling constants 
and a breakdown of the equivalence principle [31 0] . In 
this context it is particularly interesting to probe possible 
variations of the fundamental constants in the cosmolog¬ 
ical epoch of the phase transition, going from a matter- 
dominated universe to a dark energy-dominated universe, 
covering redshift ranges z = 0.5 — 5 (S3- While models 
of Big Bang nuclear synthesis probe fundamental con¬ 
stants at extremely high redshifts (z = 10 8 ) [6], the Oklo 
phenomenon (z = 0.14) [7] and laboratory atomic clock 
experiments (z = 0) [S] probe low redshifts. Absorbing 
galaxies in the line-of-sight of quasars are particularly 
suitable for investigating the range of medium-high red¬ 
shifts, for a varying fine-structure constant, a, based on 
metal absorption f9| and for a varying proton-electron 
mass ratio, p = m p /m e , based on molecular absorption 
m ■ Furthermore, unification scenarios predict that vari¬ 
ations of a and p are connected, while in most schemes 
p is a more sensitive target for varying constants m 

A variation of p may be probed through the spec¬ 
troscopy of molecules such as hydrogen (H2)j32], ammo¬ 
nia (NH3)[T3] and methanol (CH30H)[T4|. The latter 
polyatomic molecules are more sensitive to a variation in 
p under the assumption that molecular reduced masses, 
involving protons and neutrons, scale in a similar man¬ 
ner as p, which can be probed in a pure form in H2 
(TO] , Moreover, the radio absorption systems, where the 
polyatomic absorbers can be found, are more rare and 
are currently only found at the lower redshifts z = 0.68 


mm and z = 0.89 [14]. Molecular lines observed at 
z = 6.34 pT] are of insufficient spectral quality to con¬ 
strain /i-variation. Conveniently, the ultraviolet spec¬ 
trum of the H2 Lyman and Werner bands can be in¬ 
vestigated with large ground-based optical telescopes for 
absorbers at redshifts z > 2. The current sample of 
H 2 -based measurements covers a redshift interval from 2 
to 3 LIU US], with the highest redshift object Q0347—383 
at z = 3.025 EDI [2D. 

Here we constrain variations in p at substantially 
higher redshift by analyzing an H2 absorber at z = 4.22 
along the sightline towards the background (z = 4.42) 
quasar PSS J1443+2724 [22] ■ This step to higher redshift 
is challenging for several reasons. Firstly, more distant 
quasars are typically fainter, making initial discovery of 
the H2 absorption more difficult and requiring longer in¬ 
tegration times for a high-quality spectrum with which to 
constrain A p/p. Secondly, absorption lines from neutral 
hydrogen (Hi) in the intergalactic medium are more nu¬ 
merous at higher redshifts, complicating analysis of the 
H2 transitions. The Lyman-series transitions from the 
many unrelated Hi clouds in the intergalactic medium 
form a characteristic ‘forest’ of broader spectral features 
against which the H 2 have to be identified and analyzed 
(see Fig. [l]) . However, the ‘comprehensive fitting’ method 
of simultaneously treating this Hi forest and the H2 ab¬ 
sorption, developed previously and documented exten¬ 
sively ESI EDI E31 El], is employed here to reliably meet 
this challenge. 

A portion of the J1443+2724 spectrum is shown in 
Fig. [TJ To create the spectrum, we make use of an 
archival dataset obtained in 2004 [22] and a new dataset 
obtained in 2013 (Program ID 090.A-0304) which com¬ 
prise, respectively, 7.3 and 16.7 hours (5 and 12 ex¬ 
posures) observation time with the Very Large Tele¬ 
scope/the Ultraviolet and Visual Echelle Spectrograph 
(VLT/UVES). For weak light sources as the present 
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FIG. 1. Part of the final, fitted spectrum of J1443+2724. The sticks indicate the 3 velocity components of each Ht transition 
from the rotational levels J = (0 — 4). The residuals are shown above the spectrum for the regions fitted. The broad features 
that surround and overlap the fitted H 2 transitions are multiple, unrelated Hi Lyman-series absorption lines arising in other 
absorbers along the line of sight. 


TABLE I. Results of fitting models of increasing complexity 
to the spectrum. The two velocity features (see Fig. [2| re¬ 
quire fitting at least 2 velocity components (VC) but as it can 
be seen from the goodness-of-fit measure yy (where v is the 
number of degrees of freedom; here v ~ 4600) models with 
3, 4, and 5 VC replicate the data better. For each model, a 
resulting A/r//x value is shown in the last column. 


Number of VCs 

xl 

An/fj, [xl0“ 6 ] 

2 

1.396 

-8.7 ±5.2 

3 

1.161 

-6.7 ±5.4 

4 

1.139 

-8.3 ±5.4 

5 

1.133 

-2.9 ±5.2 



quasar (visual magnitude V ~ 19.4) long integration time 
is required to reach signal-to-noise ratio (SNR) of >30, 
especially if high spectral resolution is required. Wave¬ 
length calibration in the UVES instrument is achieved 
by observing a reference ThAr lamp. After each night 
of observations (in 2004) or, preferably, immediately af¬ 
ter each quasar exposure (in 2013) a spectrum of the 
ThAr lamp was recorded using identical settings to those 
of the science exposure. The ThAr spectrum was used 
by the UVES data-reduction pipeline to create pixel- 
vs.-wavelength maps and apply them to the quasar ex¬ 
posures. The wavelength scales of the spectra were 
converted to vacuum-heliocentric values. Following the 
same procedures as in [HI [23], [243 j hux-extracted individ¬ 
ual exposures were then resampled, scaled and merged 
to a final ID spectrum that extends from 474 to 792 nm, 
with the H 2 transitions detected at 484-581 nm. The 
SNR at shorter wavelengths is entirely dominated by the 
data from 2004 as the observations in 2013 were affected 
by stray light from the full moon. 

As shown in Fig. [2] each H 2 transition in the 
J1443+272 spectrum has two distinct velocity features. 
In the comprehensive fitting method employed here, 
these two features in all H 2 transitions are fitted sirnul- 


FIG. 2. (a) Composite residual spectra of 31 transitions 

shown for the 3, 4, and 5 VC models, (b) An H 2 absorp¬ 
tion profile in the J1443+2724 spectrum. The sticks indicate 
the positions of velocity components for the 3, 4, and 5 VC 
models. 


taneously (along with the many broader, unrelated Hi 
forest features). With this approach, one makes use of 
the known molecular physics of H 2 , thereby reducing the 
number of free parameters in the fit. In particular, tran¬ 
sitions from the same ground rotational level J are fitted 
using a single parameter for column density, Nj , while, 
for the same velocity component (VC), transitions from 
all J levels are tied in terms of redshift, z, and line-width, 
b. The intrinsic intensities of the H 2 absorbing transi¬ 
tions are fixed to the oscillator strengths known from the 
molecular physics database [23|. In this way, the best 
fit is achieved by simultaneously combining information 
from multiple spectral regions. In the case of J1443+272, 
we selected 60 spectral regions containing a total of 89 
H 2 transitions up to J = 4, among which 17 are from 
the Werner band [25] . This approach, along with the 
simultaneous fitting of the Hi forest lines, is the same 
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as detailed previously (e.g. [HEIM!)- The x 2 minimiza¬ 
tion is performed using a Voigt profile fitting program 
VPFIT9.5 [26]. To fit each transition i. a Voigt profile is 
created from a threefold convolution of a Lorentzian pro¬ 
file defined via a damping parameter T*, a Gaussian pro¬ 
file describing the thermal and turbulent velocities in the 
gas, and an instrumental profile. For the z-th transition 
of H 2 detected at redshift z a b s , the recorded wavelength 
is expressed as: 

K = A° (1 + Zabs)(l + -fr-i—~)j (1) 

where A? is the corresponding laboratory wavelength, 
and Ki is a coefficient that quantifies the sign and magni¬ 
tude of its sensitivity to a varying fj, [12] , and where that 
variation is parametrized by A////X = (/r z — AhauV/hab- 

The two velocity features per transition can contribute 
towards better precision of A/////. However, previous 
studies have shown that an inadequate fit of H 2 features 
can produce a spurious A/i//z value [223 ■ I n particular, 
underfitted velocity profiles are more prone to such er¬ 
rors. To avoid that, additional VCs are included in the 
profile to check for a more complex velocity structure 
than can be appreciated by eye. The additional compo¬ 
nents can be rejected by VPFIT as statistically unjustified 
or they can remain and improve the model. Table [I] con¬ 
tains the results of fitting increasingly complex models to 
the H 2 spectrum of J1443+272. Fig. [2] displays compos¬ 
ite residual spectra of different models. The composite 
residual spectra are created by combining a number of 
residuals of individual transitions, which are aligned in 
velocity/redshift space [23] . While according to the re¬ 
duced x 2 value, a model with 4 or 5 VCs initially might 
appear preferable to one with only 3 VCs, fitting the for¬ 
mer more complex models was not stable: vpfit rejected 
the additional VCs in lower-J transitions while retaining 
them in higher-J transitions, thereby departing from a 
physically plausible, self-consistent model. Fig. [2] shows 
that adding a second VC to the left feature of the absorp¬ 
tion profile reduces the composite residuals, but that the 
residuals for the right feature are barely affected when 
adding a third VC. 

In all the models, the H 2 absorption feature on the 
right requires fitting of at least one very narrow com¬ 
ponent. This weak component exhibits an unusually 
high relative column density for J = 2,3 levels, and a 
width of 0.2 km/s, corresponding to a kinetic tempera¬ 
ture of 5 K, hence lower than Tcmb = 14 K at z = 4.22. 
While this low b value results from a best-fit model and 
hence is favored statistically, we performed further test¬ 
ing by imposing 6 m ; n parameters corresponding to tem¬ 
peratures including those in the expected kinetic regime 
of 50-100 K (b ~ 0.6 —0.9 km/s). These tests, performed 
for 3 VC and 4 VC models demonstrate that resulting 
values for A/z//z do not critically depend on the narrow¬ 
ness of this velocity feature [25]. A spectrum of a higher 


TABLE II. Results of fitting a 3 VC model to the data. 

Sabs b [km/s] log N±aiogN 

from J = 0 to J = 4 

4.2237292(6) 1.74T0.06 17.04i0.06, 17.40i0.05, 15.86i0.09, 
14.81i0.04, 12.99i0.84 

4.2240440(11) 0.42i0.18 17.49i0.03, 17.87i0.02, 17.34i0.02, 
17.09i0.03, 13.80i0.33 

4.2241975(20) 0.20i0.08 15.01i0.33, 16.17i0.10, 15.77i0.09, 



FIG. 3. Statistical constraints obtained from the fidu¬ 

cial 3 VC model, from various tests performed on it, and from 
the 4 VC and 5 VC models. The green points refer to fitting 
low- and high-J transitions with the z and b parameters tied 
within each group but not between the two groups. The left 
and the right tests refer to deriving A/r//r constraints from the 
two velocity features of H 2 separately (see Fig. [2j|. The or¬ 
ange point was obtained after the spectrum obtained in 2004 
has been corrected for long-range wavelength distortions. 

SNR and higher resolution might help in explaining the 
composite residual excess and the N, b values obtained 
here. 

Based on the internal consistency of various fitting re¬ 
sults, a model with 3 VCs is selected as the fiducial model 
(see Table |II] for fitting results) and is further tested for 
possible systematics. The 4 and 5 VC fits are used to 
verify the results of the fiducial model. 

The fiducial 3 VC model delivers A/n//z = (—6.7 ± 
5-4 sta t) x 10 -6 . Fig. [5] shows how this result compares 
to constraints obtained from fitting only certain selected 
transitions, velocity components, or parts of the spec¬ 
trum. Most of the derived constraints are compatible 
with the fiducial result, demonstrating that the result 
is robust. Nevertheless, a possible discrepancy is that 
constraints from the low- and high-J transitions seem 
to move in opposite directions if fitted separately. This 
tendency is observed in the 4 and 5 VC models as well. 
If the redshift and width parameters are tied separately 
for the low-J and high-J transitions the resulting z and 
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b parameter values agree well between the two groups. 
Fitting only the small number of Werner transitions (17), 
together with their smaller spread in A) coefficients, re¬ 
sults in a large uncertainty in A/////. Fitting a spec¬ 
trum that includes only the exposures from 2004 delivers 
A [i/n = (—0.2± 6.4 sta t) x 10 -6 . This constraint, as well 
as the constraints derived from the 4 and 5 VC models 
on a full spectrum, are in agreement with the result from 
the fiducial model. Fitting a spectrum that includes only 
2013 exposures results in A/z//z = (4.0 ± 12.4 stat ) x 10 -6 . 
The uncertainty is much larger because the SNR of the 
data collected in 2013 is affected by additional noise from 
the full moon. 

The accuracy to which A/z//z can be measured strongly 
relies on accurate wavelength calibration of the quasar 
spectrum. At any wavelength in the J1443+2724 spec¬ 
trum, the error of the wavelength calibration solution is 
~ 70 m/s which translates to A/z//z of lxlO” 6 given a 
spread in A) coefficients of 0.05 ( 29] . 

An additional source of systematic uncertainty might 
be caused by wavelength-scale distortions inherent to the 
instrument. Short-range distortions with a characteristic 
shape repeating in each diffraction order of the echelle 
spectrograph were found in UVES m- We estimate its 
potential to systematically shift A/z//z value by applying 
a saw-tooth like distortion of ±100 m/s to each order 
and fitting the recombined spectrum, as in [THUS]. This 
results in a A n/n shift as small as 1 x 1(U 7 . 

Recently, indications were reported that the UVES in¬ 
strument is also susceptible to long-range wavelength dis¬ 
tortions that are able to mimic non-zero values 

at the level of several parts per million m- Their ori¬ 
gin is not clear but deviations between the light paths 
of the quasar and ThAr lamp is considered the likely 
(proximate) cause. The effect of such a miscalibration 
can be quantified by observing objects that have well- 
understood spectra, such as asteroids which reflect the 
solar spectrum and stars with solar-like spectra, known 
as solar twins P33 HUGS 122] ■ In the archive of UVES [52] . 
we found several sun-like stars that were observed in 
2004, at a similar time as the J1443±2724 quasar 834] , 
By comparing their spectra to a highly accurate Fourier- 
Transform spectrum of the sun [55] we found a correc¬ 
tion amounting to 44.9±46.8 m/s/1000 A. An extensive 
study of this particular systematic effect shows that cor¬ 
rections of this size are typical for the considered time 
period |32j . Applying the aforementioned correction on 
the wavelength scale of the spectrum obtained in 2004 
results in a A fi//i shift of —2.8 x 1(U 6 if compared to 
the uncorrected constraint. The uncertainty on this cor¬ 
rection translates into a systematic A fj/fi uncertainty of 
1.6 x 10 -6 . We consider this shift a representative cor¬ 
rection for the long-range wavelength distortions of the 
total spectrum since 84% of the fiducial A/z//z accuracy 
is gained from the 2004 exposures. 

Creating a ID spectrum from multiple exposures and 
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FIG. 4 . The non-relativistic mass density parameter and 
the dark energy density parameter Da as functions of redshift, 
with a corresponding lookback time axis on top. Results from 
various A fi/fi studies: the diamond markers refer to the mea¬ 
surements of rotational transitions of methanol and ammonia 
EMU, while A/z/p values derived from H 2 observations are 
denoted by circles [ 1844211123112411361137 ] . 


overlapping echelle orders involves redispersion onto a 
single wavelength grid. The rebinning causes correlation 
in flux and flux uncertainty values between adjacent pix¬ 
els. Varying the dispersion bin size from the default value 
of 2.5 krn/s by ±0.1 km/s results in an average A\x/\x 
shift of ±4 x 10" 6 . 

The spectra were obtained using 1.0 and 0.8" slitwidths 
in 2004 and 2013, respectively. This corresponds to re¬ 
spective spectral resolutions of R ~ 38 700 and ~ 48 500 
or Gaussian instrumental profiles with the width of 
<7inst = 3.3 and 2.6 km/s. However, the profile to be con¬ 
volved with the fitted line profile will be up to ~10 % nar¬ 
rower because the quasar light is concentrated towards 
the centre of the slit. The fiducial A/j//i constraint was 
derived from the combined 2004 and 2013 spectra using 
an instrumental profile of dinst = 2.6 km/s. Varying the 
dinst in the range from 2.6 to 2.9 km/s resulted in de¬ 
viations from the fiducial A /i/fi constraint as large as 
3 x 1CT 6 . 

In summary, the constraint on A fi/ix derived from 
this quasar spectrum shows sensitivity to these system¬ 
atic effects. All the contributions combined (added in 
quadrature) constitute the total systematic error budget 
of a sys = 5.3 x 10 -6 . After correction for the long-range 
wavelength distortions, our fiducial constraint therefore 
becomes A/j,//i = (—9.5 ± 5.4 sta t ± 5.3 sys ) x 10 -6 . 

This result constrains a possible variation of the 
proton-electron mass ratio by means of the H 2 spectro¬ 
scopic method for the highest redshift so far. Fig. ^dis¬ 
plays a comparison with previous data, mainly showing 
results derived from H 2 absorbers toward quasars in the 
redshift interval z = 2 — 3, as well as the more constrain- 
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ing results from radio astronomy of polyatomic molecules 
for z < 0.9. With the analysis of the absorber toward 
J1443+2724 the window z > 4 is opened. A further com¬ 
parison can be made with a result from radio astronomy 
of a lensed galaxy probing a (7-6) rotational transition 
in CO and a fine-structure transition in atomic carbon 
( 3 4 P 2 - 3 Pi) at 2 = 5.2 corresponding to 12.7 Gyr lookback 
time [38] : that study probes the combination of dimen¬ 
sionless constants F =1 a 2 / fi and yields A F/F < 2x 10~ 5 

E9- 

The presented H 2 constraint signifies a null result. It 
will assist in setting boundaries to various theories de¬ 
scribing physics beyond the Standard Model mm , as 
well as on the ACDM standard model of cosmology [ID] ■ 
The densities of matter (including dark matter) f2 m and 
dark energy Qa> and their ratio (see Fig. [4]) are important 
parameters in the models [3], and the present work cov¬ 
ers a wide set of values: in the interval z = 0 — 4.22 dark 
energy covers 12 a = 0.68 to 0.015. While a number of cos¬ 
mological scenarios suggest varying constants El HI] no 
quantification of the rate of change is predicted, except 
for some model-dependent scenarios involving screening 
and not holding for cosmological time scales [4T]. Thus, 
while a clear threshold for new physics is lacking, the 
observational and experimental targets are positioned to 
set ever tighter constraints on varying constants. The 
present study pushes a tight constraint on a varying to 
lookback times of 10% of the age of the universe today. 
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FITTED TRANSITIONS 


In Table [Sl| the 89 fitted H 2 transitions are listed, grouped by the ground J level. Of these transitions, 17 are from 
the Werner band. The fitted regions are displayed in Fig. [ST| to [S4] 


TABLE SI. 89 fitted H 2 transitions in 60 regions of the J1443+2724 spectrum. 


J level Number of tr. Transition names 


0 

1 


2 


3 


4 


7 W3R(0), L14R(0), L10R(0), L8R(0), 

L7R(0), L1R(0), L0R(0) 

20 L16P(1), L15P(1), L15R(1), W3R(1), 

L12R(1), L12P(1), W2Q(1), L10R(1), 
L10P(1), L9R(1), L9P(1), L8P(1), 
L4R(1), L4P(1), L3R(1), L3P(1), 
L2R(1), L2P(1), L1R(1), L1P(1) 

27 L16R(2), L15R(2), L13R(2), L12P(2), 

W2R(2), W2Q(2), W2P(2), L11R(2), 
L11P(2), L10P(2), W1R(2), W1Q(2), 
L9R(2), L9P(2), L8R(2), W0Q(2), 
L7R(2), L6P(2), L5P(2), L4R(2), 
L4P(2), L3P(2), L2R(2), L2P(2), 
L1R(2), L1P(2), LOP(2) 

26 L15R(3), W3P(3), L13R(3), L13P(3), 

L12R(3), W2R(3), L12P(3), W2Q(3), 
W2P(3), L11P(3), L10R(3), L10P(3), 
W1R(3), L9R(3), W0R(3), W0Q(3), 
L6R(3), L5P(3), L4R(3), L4P(3), 
L3R(3), L3P(3), L2R(3), L1R(3), 
L0R(3) 

9 L17P(4), L16P(4), L15R(4), L13P(4), 

L12R(4), L11P(4), W1Q(4), L9P(4), 
W0Q(4) 


CONTINUUM ERROR 

One possible source of uncertainty in A/i//r can in principle be caused by incorrect treatment of the continuum; 
we address it here. One can assume that, for 60 fitted sections (containing the 89 H 2 lines), the continuum error 
is random, with the error in the slope of the continuum over any individual H 2 section (of width ~100km/s) being 
Gaussian-distributed with a sigma of approximately 0.05 [norm, flux units] / 100 km/s = 5x 10~ 4 per km/s. Individual 
H 2 transitions are ~40km/s wide in this absorber, so this corresponds to an incorrect continuum slope of ~2xl0~ 3 
across the transition. To a reasonable approximation, this should impart a spurious shift to the fitted centroid of 
0.002 times the width, i.e. 40 km/s again, or a shift of ~ 80 m/s for an individual line. Some of the 89 lines are 
blended together and therefore share the same spurious shift. Thus, effectively a sample of ~70 transitions reduces 
the per-line error of 80m/s down to a total error budget of 80 /a/ 70 ~ lOm/s, which is well below our systematic or 
random error budgets. 

The global quasar continuum assumed initially is just a nominal starting guess which is refined by introducing 
local continua to each fitting region. These local continua are either constants or straight lines so that degeneracies 
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with broad Ly-a lines would be avoided. Therefore, errors in the global continuum level in each fitting section are 
marginalized over by fitting a local continuum with a slope and also fitting the broader HI Ly-a lines. 


VARYING THE MINIMUM LINEWIDTH PARAMETER 


The minimum lincwidth paramater & m j n in VPFIT is user-defined. Throughout the paper, we use a 6 m i n = 0.2km/s. 
In the case of J1443+2724, one of the H 2 velocity components reaches this limit. This is likely unphysical because 
&H 2 = 0.2km/s corresponds to T ~5K, while the temperature of the Cosmic Microwave Background radiation at 
2 = 4.224 corresponds to T = 14 K or &h 2 = 0.3km/s. Furthermore, the kinetic temperature in cold molecular 
clouds is typically around ~ 100 K or &H 2 = 0.9km/s. In Fig. S5 we show how imposing these values as 6 m i n affects 


Afi/fi. Overall, the resulting A/_x//i values are in agreement with the fiducial A\x/\x limit from the 3VC model with 
& m in = 0.2km/s. However, the goodness of fit parameter increases with increasing 6 m i n values. For example, a xt 
of 1.161 results from the 3VC model with 6 m ; n = 0.2km/s, and xt = 1-189 from the same model with b m - m = 0.9km/s. 
This indicates that the former model is statistically preferred over the latter. 


SUPERCALIBRATION 


Long-range wavelength distortions in UVES can be quantified by means of a so-called supercalibration method 
which is based on observations of objects with well understood spectra [19, 32]. In particular, these objects can be 
asteroids or solar twins which both exhibit solar-like spectra. Their observations include the usual ThAr calibration 
routine, used also for the quasar, and the resulting spectra are compared to a highly accurate solar spectrum obtained 
by an independent instrument, namely Fourier Transform Spectrometer (FTS). The comparison provides means to 
verify correctness of the ThAr calibration and to make adjustments if it is found to be flawed. 

In Table [S2} observational details are provided regarding the J1443+2724 observing runs both in 2004 and 2013. 
The service mode J1443+2724 observations were conducted during 2-3 and 17-19 of March, 2004. Note that most 


of the total SNR of the J1443+2724 spectrum originates from this dataset (see Table S31. Even though dedicated 
asteroid or solar twin observations were not carried out along with the quasar observations, incidentally, on the ESO 
archive [33] we could find several sun-like stars observed during that time period under a different program; their 
observational details are collected in Table |S4| It has to be noted that even though the same arm of the spectrograph 
was used, the settings of the quasar observations and those of the sun-like stars were not identical. By applyin g th e 
supercalibration method described above, we found corrections amounting to 44.9±46.8 m/s/lOOOA (see Fig. S6l. 
Here we assume that the distortions remain semistable over time periods as long as several weeks, and that the 
slopes obtained from the solar twins hold also for the main target. The resulting A/z//x correction is quoted in the 
manuscript. 

Asteroids Juno and Flora, and the solar twin HD126525 were observed along with J1443+2724 in March 2013. 
The instrumental settings of the quasar observations were maintained during the observations of calibration sources: 
a 544 or 520 nm central wavelength, 2x2 CCD binning, 0.8" slit. For this sample of solar-like spectra, a significant 
deviation from perfect calibration was found with an average slope value of the correction amounting to 297.6 ± 122.7 
m/s/lOOOA. In a study of H 2 toward B0642—5038, a miscalibration of similar size yielded a A/i//r correction of 
— 10 x 10~ 6 [24]. Fitting the uncorrected 2013 subspectrum results in = (4.0 ± 12.4) x 10 -6 for the 3VC model. 
If a correction is applied a A\if \i of (—9.3 ± 12.4) x 10 -6 is obtained. Because of a low SNR, it is difficult to reliably 
estimate the uncertainty of this shift. Hence, we do not correct the present 2013 subspectrunr of J1443+2724. 
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FIG. SI. Spectrum of the J1443+2724 quasar with the fitted H 2 absorption lines at redshift 2 : = 4.224 (part 1 of 4). The 
vertical tick marks indicate positions of the velocity components for the J = (0 — 4) transitions in the 3 VC model. 
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FIG. S2. Spectrum of the J1443+2724 quasar with the fitted H 2 absorption lines at redshift 2 = 4.224 (part 2 of 4). 
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FIG. S3. Spectrum of the J1443+2724 quasar with the fitted H 2 absorption lines at redshift 2 = 4.224 (part 3 of 4). 
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FIG. S4. Spectrum of the J1443+2724 quasar with the fitted H 2 absorption lines at redshift 2 = 4.224 (part 4 of 4). 
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FIG. S5. Results of fitting the J1443+2724 spectrum with different minimum width parameter fo m in- For the H 2 molecule, the 
Cosmic Microwave Background radiation at z = 4.224 corresponds to T = 14 K. The kinetic temperature in cold molecular 
clouds typically amounts to 50 — 100 K. 



Epoch of observation (2004) Epoch of observation (2013) 

FIG. S6. (a) Supercalibration results from 2004. The average correction value amounts to 44.9±46.8 m/s/lOOOA. (b) Super¬ 
calibration results from 2013. The average correction value amounts to 297.6±122.7 m/s/lOOOA. 




























TABLE S2. ESO archival data of the J1443+2724 observations with VLT/UVES (program IDs: 072.A-0346(B) and 090.A- 
0304(A)). The archival exposure name contains the date and time of the observations. 


Archival exposure name 

Central wavelength 
[nm] 

CCD binning 

Slitwidth 

["] 

Integ. time 

M 

UVES.2004-03-02T07:52:43.078.fits 

580 

2x2 

1.0 

5225.00 

LIVES. 2004-03-03T08:12:39.143.fits 

580 

2x2 

1.0 

5016.00 

UVES.2004-03-17T08:14:13.810.fits 

580 

2x2 

1.0 

5091.00 

LIVES. 2004-03-18T06:53:02.273.fits 

580 

2x2 

1.0 

5225.00 

UVES.2004-03-19T06:35:36.823.fits 

580 

2x2 

1.0 

5800.00 

LIVES. 2013-03-28T05:15:58.805.fits 

544 

2x2 

0.8 

4800.00 

UVES.2013-03-28T06:39:39.743.fits 

544 

2x2 

0.8 

4800.00 

UVES.2013-03-28T08:02:50.876.fits 

544 

2x2 

0.8 

4800.00 

LIVES.2013-03-29T05:22:20.681.fits 

544 

2x2 

0.8 

4800.00 

LIVES.2013-03-29T06:45:12.579.fits 

544 

2x2 

0.8 

4800.00 

LIVES.2013-03-29T08:10:27.130.fits 

544 

2x2 

0.8 

4800.00 

UVES.2013-03-30T05:03:29.905.fits 

520 

2x2 

0.8 

4800.00 

UVES.2013-03-30T06:28:32.340.fits 

520 

2x2 

0.8 

4800.00 

UVES.2013-03-30T07:51:l 1.513.fits 

520 

2x2 

0.8 

6000.00 

UVES.2013-03-31T05:03:05.640.fits 

544 

2x2 

0.8 

4800.00 

UVES.2013-03-31T06:25:58.450.fits 

544 

2x2 

0.8 

4800.00 

UVES.2013-03-31T07:49:l 1.563.fits 

544 

2x2 

0.8 

6300.00 


TABLE S3. Representative signal-to-noise (SNR) values measured in the continuum of the J1443+2724 spectrum. The 
redshifted H 2 transitions are found in the range between 484 and 581 nm. The total spectrum is a combination of two datasets 
obtained independently in 2004 and 2013. The SNR at shorter wavelengths is entirely dominated by the data from 2004 as the 
observations in 2013 were affected by the full Moon. 


A [nm] 

SNR to t 

SNR2004 

501.5 

25 

23 

528.5 

30 

26 

560.0 

48 

39 


TABLE S4. ESO archival data of several sun-like stars observed in 2004 (program IDs: 072.D-0739(A) and 072.B-0179(A)), at 
a similar time when also part of the J1443+2724 data were taken. The archival exposure name contains the date and time of 
the observations. 


Target 

Archival exposure name 

Central wavelength 
[nm] 

CCD binning 

Slitwidth 

["] 

HD140538 

UVES.2004-03-04T08:11:03.429.fits 

860 

2x2 

1.0 

HD140538 

LIVES. 2004-03-04T08:11:06.685.fits 

390 

2x2 

1.0 

Hip64345 

UVES.2004-03-12T07:14:53.472.fits 

564 

lxl 

0.3 

Hip67534 

UVES.2004-03-12T08:09:56.593.fits 

564 

lxl 

0.3 

Hip64459 

UVES.2004-03-25T07:43:09.439.fits 

564 

lxl 

0.3 

Hip64459 

UVES.2004-03-25T07:46:30.774.fits 

564 

lxl 

0.3 

Hip99224 

UVES.2004-03-25T09:12:28.662.fits 

564 

lxl 

0.3 

Hip99224 

UVES.2004-03-25T09:22:46.173.fits 

564 

lxl 

0.3 

Hipl02793 

UVES.2004-03-28T09:01:15.687.fits 

564 

lxl 

0.3 

Hipl02793 

UVES.2004-03-28T09:17:43.306. fits 

564 

lxl 

0.3 



